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Abstract 



The search for the top squark (ti) within the kinematic reach of Tevatron 
Run II is of great contemporary interest. Such a ti can explain the baryon 



O ' asymmetry of the universe provided 120 GeVjc? < mr < rrit. Moreover if 

00 ! _ ^ 

Q \ Am = mj^ - m~o is small, where Xi is the LSP, the dark matter relic density as 

obtained from the WMAP data may be explained via ti- LSP coannihilation. In 

" I 

^ ' this scenario the decay ti — cx\ is likely to occur with 100 % branching ratio 

H ; 

^ but for small Am the conventional di-jet + fix signal becomes unobservable. 

We propose a new search strategy based on the 2j+ signature accompanied 
by an isolated cluster of energy which arises from a decaying heavy particle with 
characteristic decay length. Our preliminary simulations with PYTHIA indicate 
that for 100 GeV/(? < mj^ < 130 GeV/c^ this signal may be observable while 
somewhat lager mj^ may still provide hints of new physics. 
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1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) [1] there are two scalar super- 
partners and tji, of the top quark which are the weak eigenstates. The mass eigenstates 
the lighter top squark ( ti ) and the heavier top squark (^2) are linear combinations of the 
weak eigenstates. Due to mixing effects in the top squark mass matrix in the weak basis 
driven by the top quark mass (m^) there may be a significant mass difference between ti 
and t2- In fact the former could very well be the next-to-lightest supersymmetric particle 
(NLSP), the lightest neutralino (x?) being the lightest supersymmetric particle (LSP) by 
the standard assumption in R-parity conserving MSSM. This happens in a wide region of 
the MSSM parameter space. In this scenario, henceforth referred to as the ti - NLSP sce- 
nario, the ti may be the only strongly interacting superpartner within the kinematic reach 
of Tevatron Run II experiments with a relatively large production cross-section. 

Additional interest in the light top-squark scenario stems from the observation that 
the MSSM can explain the baryon asymmetry of the universe via electroweak baryogen- 
esis (EWBG) provided 120 GeV/c^ < mj^ < rrit [2]. The search for ti is, therefore, a high 
priority program for the on going experiments at the Tevatron. 

The search for ti-NLSP at Tevatron Run I and LEP and, more recently, at Tevatron Run 
II produced negative results and lower bounds on rrij^. Most of the analyses [3, 4, 5] are 
based on the assumption that ti decays via the Flavor Changing Neutral Current (FCNC) 
induced loop decay, ti — > cx? [6] with 100 % branching ratio (BR). We also employ this 
assumption which is by and large valid if tan/3 > 7 where tan f] is the ratio of the vacuum 
expectation values for the two neutral Higgs bosons present in the MSSM [7, 8]. For lower 
values of this parameter the four body decay of the ti may be a competing channel [7, 8, 9] 

There are decay modes of the ti - NLSP other than the above two channels. They are 
the tree-level two body decay, ti txi and the three body decay, ti — > bWxi- The last 
modes are kinematically forbidden for small values of the mass difference Am = mj^ — m~o 
which is the main concern of this paper. 

The search for ti-NLSP are based on the jets plus missing Et channel[4, 5]. Some of 
the more recent works employed c-jet tagging by a lifetime based heavy flavour algorithm. 
These jets become softer if Am is small. As a result the efficiency of the kinematical cuts 
for suppressing the background as well as that of c-jet tagging decreases. This weakens the 
limit on m^^ from Tevatron. At Tevatron Run I the largest m-jr^ excluded was 122 GeV/c^ 
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for m~o = 55 GeV/c^. The most recent analysis by the DO collaboration at Run II [5] with 
c-jet tagging obtained the limit 150 GeV/(? for m~o = 65 GeV/c^ioi the most conservative 
cross-section after including the next to leading order (NLO) corrections [10]. 

On the other hand the LEP lower-bounds on m^^ are restricted mainly due to kinematics 
and are around 100 GeV/c^ [3]. However, much smaller values of Am can be probed in the 
cleaner environment of an — e" collider. 

The prospect of ti - NLSP search via this decay channel at Run II was investigated in 
[11]. It was observed that a large region of the mj^ - m-o parameter space corresponding to 
small Am is beyond the reach of Run II. For a given nij^ there is a minimum value of Am 
that can yield an observable signal. 

A modified strategy for ti-NLSP searches in the limit of small Am is important in its 
own right. The current interest in this search, however, is further strengthened by one of 
the cornerstones of the interface between particle physics and cosmology. A very attractive 
feature of the R-parity conserving MSSM is that the LSP (Xi), is a very good candidate for 
the dark matter (DM) in the universe required, e.g, by the Wilkinson Microwave Anisotropy 
probe (WMAP) data [12]. The DM relic density depends on the annihilation cross-section 
(thermally averaged) of a LSP pair. The coannihilation of the LSP with any other super- 
symmetric particle (sparticle) is another important mechanism for relic density production. 
This mechanism is, however, efficient only when the two coannihilating particles have ap- 
proximately the same mass. Thus in the small Am scenario ti - LSP coannihilation may 
indeed be an important mechanism for producing appropriate relic density [13]. 

The region of the parameter space of MSSM consistent with the DM relic density is 
severely constrained by the WMAP data. Nevertheless even in more restricted versions of 
the MSSM like the minimal supergravity model(mSUGRA) [14] one finds a narrow region of 
the parameter space where ti - LSP coannihilation is an important relic density producing 
mechanism[15]. 

The search for ti-NLSP with a small Am is, therefore, important irrespective of the 
question of EWBG. However, it is certainly worthwhile to check whether ti with mass in the 
quoted range preferred by EWBG can also produce an acceptable DM relic density. This 
was investigated in [6]. It was found that in a significant region of the allowed parameter 
space Am is indeed small ( see figure 7 of [6]). 

The results of [6] were illustrated by specific choices of other MSSM parameters. In 
particular EWBG in the MSSM requires certain CP violating (CPV) phases. In a certain 
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phase convention the relative phase(0^) between the higgsino mass parameter /_i and the 
SU(2) gaugino mass M2 is the most important one. EWBG usually requires 0.05 < c/)^ < 
1. However, various uncertainties in the calculation does not rule out a much smaller mag- 
nitude of this phase. Thus calculations by neglecting this phase seems to be a reasonable 
approximation[16] . 

It should, however, be emphasized that the signal proposed by us is fairly model inde- 
pendent and does not depend on the CPV phase or many of the other MSSM parameters 
at all as long as the BR(ti cXi ) is close to 100%. Under this assumption the size of the 
signal depends on mj^ through the production cross-section of the ti pair via the standard 
QCD processes and on m~o through the efficiency of the kinematical cuts. 

The starting point of our work is the observation that when mass difference Am is small, 
in most of the signal events, one of the c-quarks from ti pair decay is not energetic enough 
to produce a jet which may pass jet selection criteria of the experiments at Tevatron. It 
may be seen as an isolated energy deposit in the calorimeter coming from the decay of a 
heavy particle. We call it isolated cluster (IC). Thus the proposed signal consists of a c-jet 
of modest Et accompanied by missing energy and an isolated cluster. In order to reduce the 
background we require another hard jet in the signal which in most cases comes from QCD 
radiation. Our simulations show that a set of selection criteria based on the above features 
of the signal can isolate it from the SM background. 

In this work we do not consider the prospect of fully identifying the flavour of the heavy, 
isolated, decaying object because of the rather small statistics. This leads to inevitable 
backgrounds from, e.g., bb events and W/Z + jets events. However, we shall analyze at 
the generator level some important characteristics of this object which has the potential 
of reducing the SM backgrounds to a manageable level. At the same time we emphasize 
that this work is only suggestive of a new approach to ti search at the Tevatron and needs 
detailed detector simulation for a more deflnitive statement and that is beyond the scope of 
this work. 

We have used Pythia (v 6.206) [17] for generation of both signal and background events 
which includes generation of the parton level events followed by the decay of the partons 
hadronization and decay of their daughters. Generation of both signal and background 
events take into account initial state radiation (ISR) and flnal state radiation (FSR). The 
cross-sections (T^g and cXcc are very large and most of these events generated with low s are 
not relevant for our analysis. To sample the bb and cc events better for our purpose and save 
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Figure 1: The figure on left sliows the distributions of jets in signal (m^-^ = 120 GeV/c^, 
m~o = 110 GeV/c^) with initial and final state radiation ON (dashed) and OFF (solid). The 
figure on right shows the distributions of jets in signal and all the major backgrounds we 
have analysed. 

computer time we have used a cut pr > 3 GeV/ c for generation of bb and cc events , where 
p^T is defined in the CM frame of the colliding partons. We have used the toy calorimeter 
simulation followed by jet formation in Pythia (PYCELL). 

1. The calorimeter coverage is |?7| < 3. 

2. A cone algorithm with AR = Ar]'^ + = 0.5 has been used for jet finding with 
EC^n > 10 GeV and | Vjetl < 3 and jets have been ordered in Et- 

3. We consider leptons {i = e,/i) with > 5 GeV, \ri^\ < 3. The lepton should be 
isolated from jets {AR{jet, i) > 0.5). 

4. For charged particles ( e, fi and charged hadrons ), we have used their generator level 
momentum as track momentum when required. 

5. For jets containing a B ot a, D hadron we have used their decay length information for 
determining the presence of a long lived particle. 
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Figure 2: The distributions of of the signal (m^^ = 120 6*6^/0^, m~o = llOGey/c^) and 
the major backgrounds are shown after pre-selection (left) and after all selection cuts except 
the one on (right). 

A quark or a gluon from (mainly) FSR is seen as a jet and in most cases this jet appears 
to have the highest Et- This prompts us to consider a rather unusual signature for the signal 
events as mentioned below. 

The backgrounds, particularly bb and cc events have very large cross- sect ions and hence 
we need to generate a large number of events and retain only a small fraction of them which 
pass pre-selection for detailed analysis. We have used the following pre-selection criteria: 

1. Event should have only two jets: Nj^t = 2 (see figure 1) 

2. Events with isolated leptons are rejected. 

3. One of the jets should contain a long lived particle {B or, D hadron) and henceforth 
called matched-jet (MJ). 

4. Event should have an isolated cluster resulting from the decay of a i? or D hadron 
such that AR{jet, IC) > 0.5. The direction of the isolated cluster is defined to be the 
direction of the decaying B ot D hadron. In the final selection this cluster has to be 
identified as the signature of a long lived particle, the criteria for which are discussed 
in detail later. 



5 



5. We assume that a b-jet with 30 GeV < E?^^ < 50 GeV is tagged with a probabihty 
efe = 0.4 and for EjT* > 50 GeV, = 0.5 where is the single 6-jet tagging efficiency 
(i.e., the ratio of the number of tagged 6-jets and the number of taggable 6-jets). 

The pre-selection efficiency for the signal events is rather small: for m-jr^ = 
120 GeV/c^, m~o = 110 GeV/c^ (A) only 7.3% events survives pre-selection; for mj^ = 
120 GeV/c\ m~o = 105 GeV/c^ (B) and mj^ = 120 GeV/c^, m~o = 100 GeV/c^ (C) the 
rates are 9.7% and 9.6% respectively. For ti, bb, cc and V + jet events {V = W, Z) rates are 
0.03%, 0.41%, 0.19%, 1% and 2% respectively (see table 1 for details). The leading order 
cross-section aX Q = rrij, where Q is the QCD scale for tit* production is from [18]. All 
background cross-sections have been computed by CalcHEP (version 2.3.7) [19] at Q = V^- 
The largest background comes from Z+jets production. The variation of the cross-section of 
this process with the QCD scale is not very severe. 




Figure 3: The distributions of Ej- of the isolated cluster are shown for the signal {frij^ = 
120 GeV/ c"^, m~o = 110 GeV / c^) and the major backgrounds. 

For final selection (rejection) of signal (background) events we demand the following: 
1. For the ffist jet we require E^^'^ > 25 GeV and |?7jeti| < 1-5 and for the second jet 

\Vjet2\ < 1.5 (Cut 1). 
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Table 1: Efficiency table for the signal and SM backgrounds to pass various selection criteria. 
Column 3 shows the efficiency for pre-selection; columns 4-8 show the efficiencies for each cut 
combined with the effect of pre-selection ; column 9 shows the efficiency for final selection, 
(see text for details) The last column shows the expected number of events to pass all 
selection(rejection) criteria for £ = 8 fh~^. 
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Table 2: Production cross-sections for different mq^ are given. Signal events surviving [Nsig = 
C.a.e) for £ = 8 fb~^ for different values of Am are shown alongwitfi the corresponding 
significance S. The mass parameters are in units of GeV/c^. 

2. Events should have > ^OGeV (Cut 2) (see figure 2) 

3. The long lived particle in the matched-jet should have decay length > 1.5 mm. (Cut 
3) 

4. The isolated cluster should be central and have a minimum E^'- < 1-5 and 
E^'-^ > 5 GeV. It should have a decay length > 0.1 mm. (Cut 4) (see figure 3) 

5. In the signal we expect the isolated cluster and the matched- jet to be approximately 
back-to-back in the transverse plane. So, the cut A0(/C, MJ) > 85° (see figure 4 ) 
rejects background, particularly W+jets and Z-hjets events. The matched-jet is most 
likely to be the leading jet in W+jets and Z+jets events whereas it is the 2nd jet in 
the signal events. We therefore select events whose E^"^ < 40 GeV (see figure 5). 

We also partially reconstruct invariant mass of the matched-jet ( M^^ ) 

using the 

charged tracks associated with it. Similarly is reconstructed for the isolated 

cluster. The cuts M^^/ < 4.5 Gel// and M/^*; < 2GeV/c'^, reject the bb events and 
also reduce W+jets and Z+jets backgrounds (see figure 6). The combined effect of 
these cuts (Cut 5) for the signal and the dominant backgrounds is shown in column 8 
of table 1. 

In table 1 the selection efficiencies for each cut (1 - 5) includes the effects of pre- select ion. 
The eighth column shows the efficiency for the final selection. 
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Figure 4: The distributions of /S.(j){IC^ MJ) (in radians) are shown for the signal {mj^ = 
120 GeV/c^, m~o = 110 GeV/c^) and the W+jets and Z+jets backgrounds after all cuts 
except Cut 5 (see text). 

Although a very high rejection factor of (10~^) is achieved for the bb events and no event 
in the simulated sample survives, this may still be dangerous as a^^i is very large. Since the 
signal events do not have spectacular signatures, it is not possible to apply more stringent 
criterion on any of the features and still retain a good signal. It may be required to exploit the 
subtle features of the matched-jet and the isolated cluster to get rid of the bb events. There 
are a few features which may be exploited to this end although it may be experimentally 
challenging: 

• Number of charged tracks associated with the isolated cluster. 

• Upper cut on the lifetime related observables for the matched-jet, expected to be a 
c-jet and the isolated cluster, expected to come from the decay of a D hadron in the 
signal. 

• presence of a in the isolated cluster which carries a significant fraction of its pt- 

• More reliable reconstruction of the invariant masses using full detector information. 
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Figure 5: The distributions of E^^-^ of the signal {mj_^ = 12QGeV/c^, m~o = 110 GeV/ c^) 
and the W+jets and Z+jets backgrounds after all cuts except Cut 5 (see text). 

Since these observables are rather delicate, they may be simulated only using detailed 
detector modelling and estimation of the effects of criteria based on these observables is 
beyond the scope of this analysis. However it may be said with some degree of confidence 
that bb and other background events may be further suppressed by using judicious choice of 
such criteria. 

Our final results are presented in table 2. In view of the LEP limits [3] we have looked 
into signals with rrij^ > 100 GeV/c^ and we have chosen mj^ — m~o (in GeV/(?) = 10, 15 
and 20. The significance is defined as 5* = Ngig/ \jNbkg where Nsig ( Ni,kg ) is the number 
of signal ( background ) events passing the selection criteria for an integrated luminosity of 
£ = 8 fb~^. For m^^ < 130 GeV jc? we may hope for a discovery while beyond that it may 
be restricted to the level of a hint. The parameter space probed by us contains a part of the 
region interesting in the context of EWBG. 

Again, it should be emphasized that this very preliminary analysis is designed to provoke 
the experimentalists to scan the region interesting from the point of view of dark matter 
relic density and EWBG. 
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Figure 6: The figure shows the distributions of (left) and (right) reconstructed 

using charged tracks associated with the matched-jet and the isolated cluster in signal (m^^ = 
120 GeV/c"^, m~o = 110 GeV / c"^) and bb events after pre-selection (see text) . 
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